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1 Introduction 



Liquid crystals are probably the best materials for experimental and the- 
oretical studies of topological defects. Variety of defects, relatively simple 
experiments in which one can observe them, and soundness of theoretical 
models of dynamics of relevant order parameters make liquid crystals unique 
in this respect. Literature on topological defects in liquid crystals is enor- 
mous, therefore we do not attempt to review it here. Let us only point the 
books [|T|, 0, H in which one can find lucid introductions to the topic as 
well as collections of references. 

Our paper is devoted to dynamics of domain walls in uniaxial nematic 
liquid crystals in an external magnetic field. Static, planar domain walls were 
discussed for the first time in ||. We would like to approximately calculate 
director field of a curved domain wall. We use a method, called the improved 
expansion in width, whose general theoretical formulation has been given in 
[|, 0. Appropriately adapted expansion in width can also be applied to 
disclination lines 0. 

The expansion in width is based on the idea that transverse profiles of 
the curved domain wall and of a planar one differ from each other by small 
corrections which are due to curvature of the domain wall. We calculate 
these corrections perturbatively. Formally, we expand the director field in a 
parameter which gives the width of the domain wall, that is the magnetic 
coherence length £ TO in the case at hand, but actually terms in the expansion 
involve dimensionless ratios £ m /Ri, where Ri are (local) curvature radia of 
the domain wall. Therefore, our expansion is expected to provide a good 
approximation when curvature radia of the domain wall are much larger than 
the magnetic coherence length. For planar domain walls the perturbative 
solution reduces to just one term which coincides with a well-known exact 
solution. As we shall see below, the improved expansion in width is not 
quite straightforward - there are consistency conditions and rather special 
coordinate system is used - but that should be regarded as a reflection of 
nontriviality of evolution of curved domain walls. Actually, several first terms 
in the expansion can be calculated without any difficulty, and the whole 
approach looks quite promising. 

In the present paper we consider the simplest and rather elegant case of 
equal elastic constants. In order to take into account differences of values 
of the elastic constants for real liquid crystals one can use, for example, the 
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following two strategies: perturbative expansion with respect to deviations 
of the elastic constants from their mean value, or the expansion in width 
generalized to the unequal constants case. In the former approach, the equal 
constant approximate solution obtained in the present paper can be used as 
the starting point for calculating corrections. The case of unequal elastic 
constants we will discuss in a subsequent paper. 

The plan of our paper is as follows. We begin with general description of 
domain walls in uniaxial nematic liquid crystals in Section 2. Next, in Section 
3, we introduce a special coordinate system comoving with the domain wall. 
Section 4 contains description of the improved expansion in width. In Section 
5 we discuss consecutive terms in the expansion up to the second order in 
£ m . Several remarks related to our work are collected in Section 5. 



2 Domain walls in nematic liquid crystals 

In this Section we recall basic facts about domain walls in uniaxial nematic 
liquid crystals [@- We fix our notation and sketch background for the 
calculations presented in next two Sections. 

We shall parametrize the director field n(x,t) by two angles Q(x,t), 

/ sin G cos $ \ 
n = sin 6 sin $ . (1) 
y cos / 

In this way we get rid of the constraint n 2 = 1. 

We assume that the splay, twist and bend elastic constants are equal 
(Kn = K22 = K33 = K). In this case Frank — Oseen — Zocher elastic free 
energy density can be written in the form 

jr = — (d a Gd a G + sin 2 Qd a §d a <&). (2) 

Our notation is as follows: a = 1,2,3; d a = d/dx a ; x a are Cartesian 
coordinates in the usual 3-dimensional space i? 3 ; x = (x a ). In formula (2) 
we have abandoned a surface term which is irrelevant for our considerations. 

In order to have stable domain walls it is necessary to apply an external 
magnetic field H M, H. We assume that H Q is constant in space and time. 
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Without any loss in generality we may take 




Then the magnetic field contribution to free energy density of the nematic is 
given by the following formula 

T m = -]pCaHl cos 2 0. (3) 

Here Xa is the anisotropy of the magnetic susceptibility. It can be either 
positive or negative. For concreteness, we shall assume that Xa > 0. Our 
calculations can easily be repeated if Xa < 0. The ground state of the nematic 
is double degenerate: G = and O = n give minimal total free energy density 
T = J- e + J- m . It is due to this degeneracy that stable domain walls can exist. 
Dynamics of the director field is mathematically described by the equation 

where 

' d 6 xT. 



7i is the rotational viscosity of the liquid crystal, and 5/5n denotes the 
variational derivative with respect to n. Equation (4) is equivalent to the 
following equations for the and <3> angles 

= KAQ - | sin(20)9 a $9 a <I> - ^ Xa H 2 sin(20), (5) 
7i sin 2 0^ = Kd a (sin 2 0<9 a $), (6) 

where A = d a d a . 

The domain walls arise when the director field is parallel to the magnetic 
field Hq in one part of the space and anti-parallel to it in another. In between 
there is a layer - the domain wall - across which n smoothly changes its 
orientation from parallel to H to the opposite one, that is varies from 
to 7r or vice versa. The angle $ does not play important role. The Ansatz 

$ = $ o (7) 
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with constant $o trivially solves Eq.(6). Then, Eq.(5) is the only equation we 
have to solve. In the following we assume the Ansatz (7), hereby restricting 
the class of domain walls we consider. It is clear from formula (2) that 
domain walls with varying $ have higher elastic free energy than the walls 
with constant $. 

Let us recall the static planar domain wall [[I], 0. We assume that it is 
parallel to the x 1 = plane. Then 

6 = 6 (a; 1 ), $0 = const, (8) 

where 

©o Lx—oo = 0, 6 \ x i_+ +O0 = 7T. (9) 

One could also consider an "anti-domain wall" obtained by interchanging 
and 7r on the r.h.s. of boundary conditions (9). Equation (5) is now reduced 
to the following equation 

KQl = ^ Xa H 2 sm(2Q ), (10) 

where ' denotes d/dx 1 . This equation is well-known in soliton theory as the 
sine-Gordon equation, see e.g., ||. It is convenient to introduce the magnetic 
coherence length £ m , 

U= (K/ Xa Hl) l/2 . (11) 

The functions 

Qoix 1 ) = 2 arctan(exp - ~ X ° ) (12) 

with arbitrary constant Xq obey Eq.(10) as well as the boundary conditions 
(9). The planar domain walls are homogeneous in the x 1 = plane. Their 
transverse profile is parametrized by x 1 . Width of the wall is approximately 
equal to £ m , in the sense that for \x l — Xq\ 3> £ m values of 0o differ from 
or 7r by exponentially small terms. 

The planar domain wall solution of Eqs.(5), (6) contains two arbitrary 
constants: <3>o and x\. The arbitrariness of $o is d ue to the assumption that 
the elastic constants are equal. Then the free energy density T is invari- 
ant with respect to $ — > $ + const. If the elastic constants are not equal 
this invariance is lost, and in the case of planar domain walls $o can take 
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only discrete values nn/2, n = 0, 1, 2, 3. The constant Xg appears because of 
invariance of Eqs.(5), (6) with respect to the translations x 1 — > x 1 + const. 

Notice that 6 (xq) = vr/2. Hence at the director n is perpendic- 

ular to Hq. In fact, the boundary conditions (9) imply that for any domain 
wall there is a surface on which uHq = 0. Such surface is called the core of 
the domain wall. The magnetic free energy density T m has a maximum on 
the core. 

The planar domain wall (12) plays very important role in our approach. 
In a sense, it is taken as the zeroth order approximation to curved domain 
walls. The trick consists in using a special coordinate system comoving with 
the curved domain wall. Such a coordinate system encodes shape and motion 
of the domain wall regarded as a surface in the space. Internal dynamics of 
the domain wall, like details of orientation of the director inside the domain 
wall, is then calculated perturbatively in the comoving reference frame with 
the function (12) taken as the leading term. 

3 The comoving coordinates 

The first step in our construction of the perturbative solution consists in 
introducing the coordinates comoving with the domain wall. Two coordinates 
(a 1 , a 2 ) parametrize the domain wall regarded as a surface in the R 3 space, 
and one coordinate, let say £, parametrizes the direction perpendicular to 
the domain wall. For convenience of the reader we quote main definitions 
and formulas below ||. 

We introduce a smooth, closed or infinite surface S in the usual R 3 space. 
It is supposed to lie close to the domain wall. Its shape mimics the shape of 
the domain wall. In particular we may assume that S coincides with the core 
at certain time to- Points of S are given by X{a l ,t) 1 where a % {% = 1, 2) are 
two intrinsic coordinates on S, and t denotes the time. We allow for motion 
of S in the space. The vectors k = 1,2, are tangent to S at the point 
X(a\t) []. They are linearly independent, but not necessarily orthogonal 
to each other. At each point of S we also introduce a unit vector p(a\t) 
perpendicular to S, that is 

pX k = 0, f = l. 
1 We use the notation = df/da k . 
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The triad (X^,p) forms a local basis at the point X of S. Geometrically, S 
is characterized by the induced metric tensor on S 

9ik = X ti X tk , 

and by the extrinsic curvature coefficients of S 

— * 

Ku = pX >u , 

where i,k,l — 1, 2. They appear in Gauss- Weingarten formulas 

4 = K t 0+ rj,.*,, ft = -gt'KuXj. (13) 

The matrix (g lfc ) is by definition the inverse of the matrix (gki), i.e. g lk gu = 
5}, and T\ k are Christoffel symbols constructed from the metric tensor g ik . 
Two eigenvalues k\,k 2 of the matrix (Kj), where K l - = g d Kij, are called 

extrinsic curvatures of S at the point X. The main curvature radia are 
defined as Ri — l/fcj. 

The comoving coordinates (cr 1 , a 2 , ^) are introduced by the following for- 
mula 

x = X{a\t)+ip{a\t). (14) 

£ is the coordinate in the direction perpendicular to the surface S. In the 
comoving coordinates this surface has very simple equation: £ = 0. We will 
use the compact notation: (cr 1 ,^ 2 ,^) = (cr a ), where a—1, 2, 3 and a 3 = £. 
The coordinates (<j a ) are just a special case of curvilinear coordinates in 
the space R 3 . In these coordinates the metric tensor (G a p) in R 3 has the 
following components: 

G33 = 1, G 3 k = Gk3 = 0, Gik = N-girNJl, 

where 

i, k,l,r — 1, 2. Simple calculations give 

VG = V9N, 

where G = det{G a p), g = det(gik) and iV = det(N ik ). For iV we obtain the 
following formula 

N=l-tKi + ±t 2 (KiKl-K*Kl). 
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Components G Q/3 of the inverse metric tensor in R 3 have the form 

G 33 = x Q3k = Qk3 = Q) Qik = ( N -l£ g rl( N -l}k^ 

where 

We see that dependence on the transverse coordinate £ is explicit, while cr 1 , a 2 
appear through the tensors g^, K\ which characterize the surface S. 

The comoving coordinates (cr a ) have in general certain finite region of 
validity. In particular, the range of £ is given by the smallest positive £o(cr l , t) 
for which G — 0. It is clear that such £ increases with decreasing extrinsic 
curvature coefficients K\, reaching infinity for the planar domain wall, for 
which Kj = 0. We assume that the surface S (hence also the domain wall) 
is not curved too much. Then, that region is large enough, so that outside 
it there are only exponentially small tails of the domain wall which give 
negligible contributions to physical characteristics of the domain wall. 

The comoving coordinates are utilised to write Eq.(5) in a form suitable 
for calculating the curvature corrections. Let us start from the Laplacian 
AG. In the new coordinates it has the form 

y/G da a \ da? J 

The time derivative on the l.h.s. of Eq.(5) is taken under the condition that 
all x a are constant. It is convenient to use time derivative taken at constant 
a a . The two derivatives are related by the formula 

d . d da 13 d 

where 

^U- = -pX, = -{N- l )\g kr X,{k + ££), 

the dots stand for d/dt\ a i. Let us also introduce the dimensionless coordinate 
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Now we can write equation (5) transformed to the comoving coordinates 
(a 1 , s) (with the Ansatz (7) taken into account): 

Equation (15) is the starting point for construction of the expansion in width. 



4 The improved expansion in width 

We seek domain wall solutions of Eq.(15) in the form of expansion with 
respect to £ m , that is 

e = e + £ m e 1 + &e 2 + ... . (ie) 

Inserting formula (16) in Eq.(15) and keeping only terms of the lowest order 
(~ we obtain the following equation 

<9 2 n 1 



()s , =^sm(2& ), (17) 

which essentially coincides with Eq.(lO) after the rescaling x 1 = £ m s. Its 
solutions 

so (s) = 2 arctan(exp(s — s )), 

essentially have the same form as the planar domain walls (12), but now s 
gives the distance from the surface S. This surface will be determined later. 
In the remaining part of the paper we shall consider curvature corrections to 
the simplest solution 

o (s) = 2 arctan(exp s). (18) 

Because already Go interpolates between the ground state solutions 0, n, the 
corrections 0^, k > 1, should vanish in the limits s — > ±oo. 

Equations for the corrections 6fc, k > 1, are obtained by expanding both 
sides of Eq.(15) and equating terms proportional to £ T fe „. These equations can 
be written in the form 

LO k = f k , (19) 
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with the operator L 

i = |l_ cos(2 e ) = |L + -L--i. (20) 

The last equality in (20) can be obtained, e.g., from Eq.(17): inserting O 
given by formula (18) on the l.h.s. of Eq.(17) we find that sin(2@o) = 
— 2 sinh s/ cosh 2 s, and cos(20 o ) = 1 — 2/ cosh 2 s. The expressions fk on 
the r.h.s. of Eqs.(19) depend on the lower order contributions 6/, / < k. 
Straightforward calculations give 

fl = d s e (K r r -^pjc), (21) 

h = - sin(20 o )0 2 + sd s Q Q K)Kl + dsQ^Kl - ^pk), (22) 

h = %(d&i ~ g^xja^) - 2sin(20 o )0 1 2 - § cos(20 o )03 
+sd s e 1 K)Kl - ls 2 d s Q K r r [{Klf - 3K}Kj) 

-^djiy/gg^dkQ,) + d s Q 2 (K; - fpj(), (23) 

and 

h = i(^0 2 - sg^X^Oi) ~ fg 3k X tk X(d,Q 2 + sKp t Qi) 

- sin(20 o ) (0| + 20!0 3 - |Of) - 2 cos(20 o )0 2 2 
+sd s Q 2 K)K{ + s 3 d s Q {{KlY + \{KlK s r f - %K^K)K\) 

-^ds&.K; ({Klf - ZK)K>) - -L dj (^ggi k d k e 2 ) 
-f^i^gK^d^) + sg^djK^d^ + d s Q 3 (K r r - %pX), (24) 

where d t = d/dt, di = d/da 1 . We have taken into account the fact that O 
does not depend on a 1 . 

Notice that all Eqs.(19) for are linear. The only nonlinear equation 
in our perturbative scheme is the zeroth order equation (17). 

It is very important to observe that operator L has a zero-mode, that is 
a function ipo(s) which quickly vanishes in the limits s — > ±oo, and which 
obeys the equation 

U = 0. 
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Inserting So (s) in Eq.(17), differentiating that equation with respect to so 
and putting s = we obtain as an identity that Lip = where 

Ms) = (25) 
cosh s 

The presence of this zero-mode is related to the invariance of Eq.(17) with 
respect to translations in s, therefore it is often called the translational zero- 
mode. Let us multiply both sides of Eqs.(19) by ipo(s) and integrate over s. 
Integration by parts gives 

/oo A /*oo A 

dsipoLQk = \ dsQkLipo = 0. 
-OO J —CO 

Hence, we obtain the consistency (or integrability) conditions 

/oo 
dsMs)fk(s) = 0, (26) 
-oo 

where fk are given by formulas of the type (21) - (24). We shall see in the 
next Section that these conditions play very important role in determining 
the curved domain wall solutions. 

Using standard methods || one can obtain the following formulas for 
vanishing in the limits s — > ±oo solutions 0^ of Eqs.(19): 

&k = G[f k }+C k (a\t)Ms), (27) 

where 

rs rs 

G[fk] = -ipo( s ) / dxipi(x)f k (x) + V»i(s) / dxip ( x )fk(x). (28) 



Here ipo(s) is the zero- mode (25) and 

^i(s) = ^(sinh S + -^-) (29) 

is the other solution of the homogeneous equation 

Lip = 0. 

The second term on the r.h.s. of formula (27) obeys the homogeneous equa- 
tion L0fc = 0. It vanishes when s — > ±oo. 



11 



The solutions (27) contain as yet arbitrary functions Ck(cr\t). Also 
X(a\ t) giving the comoving surface S has not been specified. It turns out 
that conditions (26) are so restrictive that they essentially fix those func- 
tions. The extrinsic curvature coefficients K\ and the metric will follow 
from X(a\ t). 

One can worry that G[fk], k > 1, given by formula (28) do not vanish 
when s — > ±oo because the second term on the r.h.s. of formula (28) is 
proportional to which exponentially increases in the limits s — > ±oo. 
However, the integrals 

/ dxipof k 

vanish in that limit due to the consistency conditions (26). Moreover, quali- 
tative analysis of Eq.(15) shows that fk ~ (polynomial in s) x exp(— \s\) for 
large \s\, hence those integrals behave like (polynomial in s) x exp(— 2|s|) for 
large \s\ ensuring that all G[fk] exponentially vanish when \s\ — > oo. 



5 The approximate domain wall solutions 

In this Section we discuss the approximate solutions obtained with the help 
of the perturbative scheme we have just described. We present formulas for 
0i and ©2, an equation for X(a l ,t) which determines motion of the surface 
S, as well as equations for the functions C\, C*2. 

The zeroth order solution is already known, see formula (18). This allows 
us to discuss the consistency condition with k — 1. Substituting fi from 
formula (21) and noticing that 

= — t— = ipo(s) 
coshs 

we find that that condition is equivalent to 

^pjl = K r r . (30) 

This condition is in fact the equation for X. It is of the same type as Allen- 
Calm equation [|10|| , but in our approach it governs motion of the auxiliary 
surface S. 
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Let us now turn to the perturbative corrections. After taking into account 
Allen-Cahn equation (30) we have f\ = 0. Therefore, the total first order 
contribution has the form 

e, = ^M. (3i) 

cosh s 

The second order contribution O2 is calculated from formula (28) with f 2 
given by formula (22). Using the results (30), (31) we obtain the following 
formula 

6 2 = Us)Cl{a\t)+Us)K)Kl + (32) 



where 



sinhs 



2 cosh 2 s ' 

1 s 

i J 3{ s ) — ~ s cos h s ; ^lis) m (2 cosh s) 

2 2 cosh s 

s 2 sinh s 1 f s x 2 
H -s : — / ax- 



4 cosh 2 s 4 cosh s Jo cosh 2 x 
The integral in ^(s) can easily be evaluated numerically. Due to the consis- 
tency conditions, the functions Ci, C 2 in formulas (31), (32) are not arbitrary, 
see below. 

The consistency condition (26) with k = 2 does not give any restrictions 
- it can be reduced to the identity = 0. More interesting is the next 
condition, that is the one with k = 3. Inserting formula (23) for fs and 
calculating necessary integrals over s we find that it can be written in the 
form of the following inhomogeneous equation for Ci(a\t) 

£($d - gto-XrkdtCy) - ^djiy/g^dkd) - KijKjCi 

= fX r ((^) 2 - 3*5*0 • ( 33 ) 

We have also used Allen-Cahn equation (30). Equation (33) determines C 1 
provided that we fix initial data for it. Similarly, the consistency condition 
coming from the fourth order (k = 4) is equivalent to the following homoge- 
neous equation for C2 

^(d t C 2 -g kr X, r Jcd k C 2 ) 
-jjdjiy/ggi"^) - K)K\C 2 = 0. (34) 
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The formulas (16), (18), (31) and (32) give a whole family of domain 
walls. To obtain one concrete domain wall solution we have to choose initial 
position of the auxiliary surface S. Its positions at later times are determined 
from Allen-Calm equation (30). We also have to fix initial values of the 
functions Ci,C2, and to find the corresponding solutions of Eqs.(33), (34). 
Notice that we are not allowed to choose the initial profile of the domain wall 
because the dependence on the transverse coordinate s is explicitly given. 
It is known from formulas (18), (31) and (32). Any choice of the initial 
data gives an approximate domain wall solution. Of course such a choice 
should not lead to large perturbative corrections at least in certain finite 
time interval. Therefore one should require that at the initial time £ m Ci <C 
l>£mC2 <C l,£ m Kj <C 1. The domain wall is located close to the surface S 
because for large \s\ the perturbative contributions vanish and the leading 
term 2 arctan(e s ) is close to one of the vacuum values 0, ir. 

Let us remark that Eqs.(30), (33) and (34) imply that a planar domain 
wall (Kj = 0) can not move, in contradistinction with relativistic domain 
walls for which uniform, inertial motions are possible. 

In the presented approach we describe evolution of the domain wall in 
terms of the surface S and of the functions C\,C2- These functions can be 
regarded as fields defined on S. In some cases Eqs. (30), (33), (34) can 
be solved analytically, one can also use numerical methods. Anyway, these 
equations are much simpler than the initial Eq.(5). 

The presented formalism is invariant with respect to changes of coordi- 
nates a 1 , a 2 on S. In particular, in a vicinity of any point X of S we can 

— * 

choose the coordinates in such a way that = 5^ at X. In these coordinates 
Eq.(30) has the form 

K V = Rl + iV (35) 

— * 

where v is the velocity in the direction p perpendicular to S at the point X, 
and Ri, R2 are the main curvature radia of S at that point. 

As an example, let us consider cylindrical and spherical domain walls. If 
S is a straight cylinder of radius R then R 1 = 00, R 2 = —R(t), v — R and 
Eq.(35) gives 



where _R is the initial radius. The origin of the Cartesian coordinate frame 




(36) 
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is located on the symmetry axis of the cylinder S (which is the z-axis), p 
is the outward normal to S, and s = (\/r 2 — z 2 — R(t))/£ m , where r is the 
radial coordinate in R 3 . As a 1 , a 2 we take the usual cylindrical coordinates 
z,<j). Equations (33), (34) reduce to 



7i 
K 



a , Cl _(^ 1 + ^ Cl )_^ Cl ^_L (37) 

|a«c 2 - (a;c 2 + ±dlc 2 ) - = o. (38) 

If Ci,C 2 at the initial time t have just constant values Ci(0),C 2 (0) on 
the cylinder, then 

Ci(0 = In(iib/i2(0) + ^y^i(0)> ^(*) = ]^(°)- (39) 

General solutions of Eqs.(37), (38) can be found by splitting C±, C 2 into 
Fourier modes, but we shall not present them here. 

The case of spherical domain wall is quite similar. Now S is a sphere of 
radius R and Ri = R 2 = —R, v — R. Equation (35) gives 



R(t) = jR 2 -—(t-t ), (40) 
V 7i 

Now the origin is located at the center of the sphere, s = (r — R(t))/£ m , and 
p is the outward normal to S. As o k we take the usual spherical coordinates. 
Then, Eqs.(33), (34) can be written in the form 

^ - h Gsb*'*"*^ + ™k a * c - l Ci = (41) 



and 



¥ c > ~ h (sb*^"*** + artf**) ~ h a > = °- (42) 



General solution of these equations can be obtained by expanding C±, C 2 into 
spherical harmonics. In the particular case when Ci,C 2 are constant on the 
sphere S the solutions Ck(t) have the same form (39) as in the previous case 
except that now R(t) is given by formula (40). 
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In the both cases our approximate formulas are expected to be meaningful 
as long as R(i)/£ m ^> 1. 

Because we know the transverse profile of the domain wall, we can express 
the total free energy F by geometric characteristics of the domain wall. One 
should insert our approximate solution for in formulas (2) and (3) for T e 
and jF m , respectively, and to perform integration over s. The volume element 
d 3 x is taken in the form 

d 3 x = ^ m \fGd 2 ads. 
For simplicity, let us consider curved domain walls for which 

d = = C 2 
at the time to- Straightforward calculation gives 

F- KV + 2K \?\ 



- \KU [ d 2 a^(-^ + ^ - -^V) + terms of the order ( 43 ) 

J iti -K'2 tx\t\2 

where l^l denotes the area of the surface S, and V is the total volume of 
the liquid crystalline sample. The first term on the r.h.s. of this formula is 
just a bulk term which appears because the smallest value of the magnetic 
free energy density has been chosen to be equal to —K/{2^ n ). The proper 
domain wall contribution starts from the second term. This term gives the 
main contribution of the domain wall to F. One can think about the core- 
sponding constant free energy 2K/^ m per unit area. The third term on the 
r.h.s. of formula (43) represents the first perturbative correction. It is of the 
order (£ m /-Rj) 2 when compared with the main term, and within the region 
of validity of our perturbative scheme it is small. One can easily show that 
this term is negative or zero. Hence, it slightly diminishes the total free en- 
ergy. In this sense, the domain walls have negative rigidity — bending them 
without stretching (i.e., with kept constant) diminishes the free energy. 



6 Remarks 

We would like to add several remarks about the expansion in width and the 
approximate domain wall solutions it yields. 
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1. In the presented approach dynamics of the curved domain wall in the 
three dimensional space is described in terms of the comoving surface S and 
of the functions Cf., k < 1, defined on S. The profile of the domain wall has 
been explicitly expressed by these functions, by the transverse coordinate £, 
and by the geometric characteristics of S. The surface S and the functions 
Ck obey equations (30), (33), (34) which do not contain £. In particular cases 
these equations can be solved analytically, and in general one can look for 
numerical solutions. Such numerical analysis is much simpler than it would 
be in the case of the initial equation (5) for the angle 0, precisely because 
one independent variable has been eliminated. 

2. We have used £ m as a formal expansion parameter. This may seem 
unsatisfactory because it is a dimensionful quantity, hence it is hard to say 
whether its value is small or large. What really matters is smallness of the 
corrections £ m ©i, £m®2- This is the case if £ m Ci <C 1, £^C 2 <C 1 and 
£, m Kj <S 1, as it follows from formulas (31) and (32). 

3. Notice that an assumption that S coincides with the core for all times in 
general would not be compatible with the expansion in width. If we assume 
that C\ = = C*2 at certain initial time t Q , Eq.(33) implies that C\ ^ at 
later times (unless the r.h.s. of it happens to vanish). Then, it follows from 
formulas (16), (18) and (31) that ^ tc/2 at s — 0, that is on S. 

4. In the present work we have neglected effects which could come from per- 
turbations of the exponentially small tails of the domain wall. For example, 
consider a domain wall in the form of infinite straight cylinder flattened from 
two opposite sides. Its front and rear flat sides have zero curvatures, and 
according to Eq.(35) they do not move. In our approximations the domain 
wall shrinks from the sides where the mean curvature I/R1 + I/R2 does not 
vanish. Now, in reality the front and rear parts interact with each other. This 
interaction is exponentially small only if the two flat parts are far from each 
other. We have neglected it altogether assuming the 2 arctan(e s ) asymptotics 
at large s. In this sense, our approximate solution takes into account only 
the effects of curvature. 

5. Finally, let us mention that the dynamics of domain walls in nematic liquid 
crystals can also be investigated with the help of another approximation 



scheme called the polynomial approximation. In the first paper [|TT] it has 
been applied to a cylindrical domain wall, and in the second one to a planar 
soliton. Comparing the two approaches, the polynomial approximation is 
much cruder than the expansion in width. It also contains more arbitrariness, 



17 



e. g., in choosing concrete form of boundary conditions at \s\ — > oo. On 
the other hand, that method is much simpler. It can be useful for rough 
estimates. 
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